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We  study  the  few-electron  eigenspectrum  of  a  carbon-nanotube  double  quantum  dot  with  spin- 
orbit  coupling.  Exact  calculation  are  combined  with  a  simple  model  to  provide  an  intuitive  and 
accurate  description  of  single-particle  and  interaction  effects.  For  symmetric  dots  and  weak  magnetic 
fields,  the  two-electron  ground  state  is  antisymmetric  in  the  spin- valley  degree  of  freedom  and  is 
not  a  pure  spin-singlet  state.  When  double  occupation  of  one  dot  is  favored  by  increasing  the 
detuning  between  the  dots,  the  Coulomb  interaction  causes  strong  correlation  effects  realized  by 
higher  orbital-level  mixing.  Changes  in  the  double-dot  configuration  affect  the  relative  strength  of 
the  electron-electron  interactions  and  can  lead  to  different  ground  state  transitions.  In  particular, 
they  can  favor  a  ferromagnetic  ground  state  both  in  spin  and  valley  degrees  of  freedom.  The  strong 
suppression  of  the  energy  gap  to  a  ferromagnetic  state  can  cause  the  disappearance  of  the  Pauli 
blockade  in  transport  experiments  and  thereby  also  limit  the  stability  of  spin-qubits  in  double  dots 
based  quantum  information  proposal.  Our  analysis  is  generalized  to  an  array  of  coupled  dots  which 
is  expected  to  exhibit  rich  many-body  behavior. 

PACS  numbers:  73.63.Fg,73.23.-b,73.22.-f 

I.  INTRODUCTION 

Experiments  on  few-electron  double  quantum  dots  al¬ 
low  the  measurement  and  manipulation  of  the  spin  de¬ 
gree  of  freedom  of  the  confined  electron^  Such  con¬ 
trol  is  at  the  heart  of  semiconductor-based  spintronics2-^ 
and  quantum-information  proposal^4  54  Recently,  sub¬ 
stantial  experimental  efforts  have  been  focussed  on  con¬ 
trolling  electrons  in  carbon  nanotube  double  quantum 
dots,  and  many  of  the  capabilities  previously  achieved  in 
GaAs  double  dote1  are  starting  to  be  reproduced®^1: 

These  include  the  ability  to  start  from  an  empty  dou¬ 
ble  dot  and  systematically  fill  it  with  electrons.  Since 
12  C  has  no  nuclear  spin,  carbon  based  nanostructures 
are  expected  to  reduce  hyperfine  induced  decoherence  as 
compared  with  GaAi1^. 

Furthermore,  carbon  based  materials  exhibit  richer 
physics  than  GaAs  semiconducting  materials  because  of 
the  additional  valley  degree  of  freedom131®.  In  princi¬ 
ple,  the  spin  and  valley  degrees  of  freedom  could  lead  to 
a  SU(4)  symmetry  at  zero  magnetic  field  instead  of  the 
standard  SU(2)  symmetry  in  conventional  semiconduc¬ 
tors  (see  e.g.  Ref.  US  and  references  therein).  However, 
it  has  been  recently  demonstrated1®  that  the  enhance¬ 
ment  of  the  spin-orbit  splitting  in  small-radius  nanotubes 
breaks  the  four-fold  degeneracy  of  the  single-electron 
ground  state  into  a  two-fold  degeneracy. 

In  this  work,  we  study  how  spin-orbit  coupling  and 
electron-electron  interaction  effects  are  manifested  in  the 
two-electron  spectrum  and  transport  properties  of  a  car¬ 
bon  nanotube  double  dot.  This  represents  an  extension 
of  previous  studies  on  few-electron  physics  in  a  single 
carbon  nanotube  dotP^l.  We  find  that,  despite  of  spin- 
orbit  coupling  and  the  existence  of  an  additional  val¬ 
ley  degree  of  freedom,  the  two-electron  eigenstates  can 
be  separated  in  an  orbital  part  and  a  spin-valley  part 


that  are,  to  a  very  good  approximation,  independent  of 
each  other.  The  spin-valley  part  can  be  grouped  in  ten 
antisymmetric  and  six  symmetric  spin-valley  eigenstates 
which  we  refer  to  as  multiplets. 

The  separation  of  spin-valley  degrees  of  freedom  sig¬ 
nificantly  simplifies  the  description  of  the  systems  and 
allows  us  to  draw  analogies  with  standard  GaAs  double 
dots.  Our  main  results  can  be  summarized  as  follows: 
(a)  For  dots  at  zero  magnetic  field  and  no  detuning,  each 
dot  is  populated  by  a  single  electron  and  tunneling  is  sup¬ 
pressed  because  of  Coulomb  interactions.  Thus,  interdot 
coupling  only  occurs  virtually  via  superexchange  interac¬ 
tions  that  determine  the  ground  state  symmetry.  In  this 
regime,  we  find  a  spin-valley  antisymmetric  ground  state, 
that  does  not  have  a  well  defined-spin  due  to  the  spin- 
orbit  coupling,  (b)  For  large  detuning  between  the  dots, 
double  occupation  of  the  same  dot  becomes  favorable. 
In  this  regime,  Coulomb  interactions  can  mix  higher  or¬ 
bitals  in  the  two-electron  ground  statd4®®1.  This  ad¬ 
mixture  significantly  reduces  the  energy  spacing  between 
the  multiplets  (which,  for  weakly  interacting  electrons,  is 
determined  by  the  orbital  level  spacing).  The  interplay 
between  spin-orbit  coupling  and  interaction  then  leads  to 
a  ferromagnetic  ground  state  above  a  small  critical  mag¬ 
netic  field,  since  the  Zeeman  terms  overcome  the  strongly 
suppressed  splitting  between  effective  singlet  an  triplet, 
c)  The  reduction  of  the  energy  gap  between  orbitally  sym¬ 
metric  and  antisymmetric  states  caused  by  the  Coulomb 
interaction  affects  transport  properties  through  the  dot 
and  can  lead  to  the  disappearance  of  the  so-called  Pauli- 
blockade  (suppression  of  current  through  the  double  dot 
due  to  the  Pauli  exclusion  principle)  and  might  explain 
the  absence  of  Pauli  blockade  reported  in  Ref.  HU  The 
absence  of  spin-blockade  might  affect  the  performance  of 
quantum  information  proposal  which  use  spin-qubits  in 
double  dots,  since  gate  operation  in  those  proposals  is 
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based  on  the  spin-blockade  mechanism.  The  disappear¬ 
ance  of  Pauli  blockade  might  be  prevented  by  reducing 
Coulomb  correlations  by  either  working  with  short  dots, 
or  by  covering  the  nanotube  by  large  dielectric^. 

This  paper  is  organized  as  follows.  In  the  section  II,  we 
introduce  the  microscopic  model  for  the  double  dot  and 
analyze  the  non-interacting  predictions  taking  into  ac¬ 
count  a  magnetic  field  parallel  to  the  nanotube  axis,  spin- 
orbit  couplings  and  detuning  between  the  dots.  Then,  we 
construct  a  simple  two-electron  model  that  captures  the 
interaction  effects.  This  model  is  then  compared  with 
solutions  of  an  exact  many-band  Hamiltonian  using  lo¬ 
calized  single-particle  orbitals.  In  Section  m  we  discuss 
the  energy  spectrum  of  two  interacting  electrons  in  a  dou¬ 
ble  dot  in  three  different  detuning  regimes  corresponding 
to  (i)  a  symmetric  double  dot,  with  one  electron  in  each 
dot,  (ii)  strong  detuning,  with  both  electrons  in  the  same 
dot,  and  (iii)  at  the  crossover  between  both  regimes.  We 
then  analyze  the  transport  properties  of  the  double  dot. 
Finally,  we  discuss  how  to  extend  our  low  energy  behavior 
analysis  to  serially  coupled  quantum  dots.  In  section  IV, 
we  present  the  conclusion.  Technical  details  on  the  cal¬ 
culation  of  Coulomb  matrix  elements  and  the  derivation 
of  the  rate  equations  used  for  transport  are  presented  in 
Appendices  [A]  and  [5] 
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FIG.  1:  (Color  online)  Schematic  representations  of  a  double 
quantum  dot  in  carbon  nanotube.  Top:  Blue  regions  corre¬ 
spond  to  the  potential  barriers.  Bottom:  Double  dot  potential 
at  finite  detuning  A.  Dashed  lines  schematically  represent  the 
single-particle  orbitals. 


II.  MODEL 

Single-wall  carbon  nanotubes  are  formed  by  a  single 
layer  of  graphite  called  graphene  rolled  up  into  a  cylinder. 
Graphene  has  a  honeycomb  lattice  formed  by  covalently 
bonded  carbon  atoms.  Its  electronic  properties  are  de¬ 
termined  by  the  pz  orbital  of  the  carbon  atom.  The  low 
energy  spectrum  of  graphene  consists  of  two  Dirac  cones 
located  at  the  K  and  K'  =  —K  points  of  the  graphene’s 
Brillouin  zone,  where  the  valence  and  conduction  bands 
touch.  To  characterize  the  two  Dirac  cones,  we  introduce 
the  valley  index  r  =  ±1,  where  r  =  1  corresponds  to  the 
K'  point  and  r  =  —  1  to  the  K  point.  The  behavior  of 
graphene  in  the  presence  of  an  external  potential  can  be 
described  by  an  effective  mass  approximation,  or  k  ■  p 
theory  21.  In  this  approximation,  the  envelope  wave  func¬ 
tion  for  the  A  and  B  sites  of  the  two-atom  unit  cell  in 
a  honeycomb  lattice  follows  an  effective  Dirac  equation. 
In  a  carbon  nanotube,  the  cylindrical  structure  imposes 
a  quantization  condition  that  leads  to  either  metallic  or 
semiconducting  nanotubes,  depending  on  the  orientation 
of  the  underlying  lattice  with  respect  to  the  symmetry 
axis  of  the  tube.  Here,  we  will  focus  on  the  behavior  of 
semiconducting  nanotubes. 


A.  Single  particle  spectrum  and  interactions 

In  this  section  we  follow  previous  worli12^^  derive 
the  localized  eigenstates  of  a  semiconducting  nanotube 
with  an  additional  confinement  potential  along  the  tube, 
which  is  controlled  by  external  gates.  We  describe  the 
confinement  potential  of  each  dot  by  a  square  well12  22 
(see  Fig.  [l]).  The  form  of  the  confinement  potential  will 
not  affect  our  results  qualitatively  and  we  note  that  for 
a  single  dot,  the  results  for  parabolic  confinement  and 
square  well  are  in  good  agreement18  20.  The  external 
potential  leads  to  a  discrete  set  of  bound  states.  Taking 
(  as  the  direction  of  the  nanotube  axis  and  ip  as  the 
angle  perpendicular  to  the  nanotube  axis,  we  can  write 
the  single  particle  Hamiltonian  as 

H0  =  -ihv(Ta1^dip  +  a2d<:)  +  V(C),  (1) 

where  v  is  the  Fermi  velocity,  cq  are  the  Pauli  matri¬ 
ces  operating  over  the  sublattice  space,  and  H(C)  is  the 
external  potential  that  describes  one  or  two  dots.  The 
eigenstates  are  determined  by  matching  the  solutions  for 
the  dot  and  barrier  regions,  which  are  of  the  form: 

4'r'K’fc(</5,  C)  =  ei(-KRv+kC)  ^  Z^k  ^  .  (2) 

Here  n,  k  denote  the  wave  vectors  around  and  along  the 
tube,  zj,  k  =  ±(tk  —  ik)/Vn2  +  k2,  and  the  energy  is 
given  by  EK  k  =  ±Hv\/k2  +  k2.  Solving  the  effective 
Dirac  equation  for  Vid(C)j  which  is  0  for  0  <  (  <  L 
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and  Vg  otherwise,  leads  to  a  quantization  condition  for 
the  longitudinal  momentum  modes  kn  of  the  localized 
states,  where  n  denotes  the  band  inde:xP^. 

So  far,  the  K  and  K'  solutions  are  degenerate  and 
independent  of  spin,  er  =//,  leading  to  a  four-fold  sym¬ 
metry.  However,  this  symmetry  is  broken  by  spin-orbit 
coupling  corrections1'  23  24  and  by  a  constant  magnetic 
field  B  along  the  nanotube  axis  £.  The  spin-orbit  cou¬ 
pling  and  the  presence  of  a  magnetic  field  modify  the 
quantization  condition  in  the  (p  direction, 

k  =  k0  +  $ab/$0R-  sASo/(hv)  (3) 

kq  =  t/3R.  (4) 

Here,  s  =  ±1/2  is  the  quantum  number  correspond¬ 
ing  to  the  spin  operator  parallel  to  the  nanotube 
(5c|t>  =  1/2  It)  and  |J.)  =  -1/2 1+));  Aso  « 
ImeV /i?[nm]  is  the  energy  splitting  due  to  spin-orbit  cou¬ 
pling,  $ab  =  BnR2  is  the  Aharonov-Bohm  flux  through 
the  nanotubc25!  and  4>o  =  hc/\e\.  The  magnetic  field 
also  leads  to  a  spin  Zeeman  term  Hz  =  shtu ,  where 
u)  =  \e\gB/(2moc)  is  the  Zeeman  frequency  in  terms  of 
the  gyromagnetic  constant  g,  the  electron  mass  mo,  and 
the  speed  of  light  c.  The  sign  convention  for  the  spin 
used  here  is  opposite  to  the  one  of  Ref.  (TSJ  Note  that  in 
Eq.  @  we  only  consider  the  lowest  mode  in  the  trans¬ 
verse  direction,  since  excitations  involve  energies  of  about 
hv/R  which  are  much  larger  than  longitudinal-excitation 
energies  or  Coulomb-interaction  effects  as  long  as  R  <C  L. 

An  example  of  the  singe  particle  energy  spectrum  of  a 
single  dot  is  shown  in  Fig.  [2]  Note  that  we  measure  the 
energy  with  respect  to  the  center  of  the  gap,  so  that  the 
dominant  part  of  the  single-particle  energy  is  constant 
and  given  by  hvn  ss  220meV/i?[nm].  Generally,  disorder 
or  the  confinement  potential  itself  can  lead  to  intervalley 
coupling.  However,  for  a  noticeable  effect,  the  potential 
must  change  on  the  scale  of  the  nearest  neighbor  lattice 
spacing  between  carbon  atoms  oo  =  1.4  A.  The  experi¬ 
ment  of  Ref.  m  shows  only  a  tiny  valley  mixing,  and  we 
neglect  intervalley  scattering  in  this  work. 

The  longitudinal  wave  vector  depends  indirectly  on 
spin-valley  quantum  numbers  and  on  magnetic  field  and 
these  effects  are  included  in  the  multi  band  calculations. 
However,  since  k0  $ab/($oR),  Aso/(2?iw),  kn,  the 
single  particle  energies  can  be  significantly  simplified, 
yielding 

En,a  «  E°{Vg)  +  Ea  (5) 

Ea  =  B{t gorb  ±  2gspins)  Asots  .  (6) 

Here  we  have  introduced  a  single  quantum  number  a  = 
(s,r)  to  describe  the  spin  valley  degrees  of  freedom 
(a  =f  K,  t  K,l  K',l  K').  The  confinement  energy 
En(Vg)  =  hvy/i^  +  k^  is  the  single-particle  spectrum  of 
a  dot  with  potential  depth  Vg  ignoring  magnetic  field, 
spin,  and  spin-orbit  dependences;  /ispm  =  hui/(2B)  = 
?i|e|<7/(4TOoc);  and  gori,  =  (iTR\e\/hc).  Within  this  ap¬ 
proximation  the  orbital  part  of  the  single-particle  states 
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FIG.  2:  (Color  online)  Single  particle  spectrum  as  a  function 
of  B  for  L  =  70  nm,  Vg  =  78  meV,  and  R  =  2.5  nm  and 
Aso  <  0.  Solid  curves  correspond  to  EnK't,  dashed  curves 
correspond  to  E„k±,  dash-dotted  curves  correspond  to  EnKi^. 
and  dotted  curves  correspond  to  EnK/±,  where  n  =  1,2,3  is 
the  band  label. 


separate  from  the  spin-valley  part.  Figure  [2]  illustrates 
that,  for  the  parameters  studied  in  this  work,  spin-valley 
splitting  is  basically  the  same  for  all  longitudinal  bands, 
and  Eq.  <[5j)  provides  a  good  description  of  the  single¬ 
particle  spectrum  of  a  single  dot. 

Next,  we  consider  a  biased  double-dot  system, 
schematically  presented  in  Fig.  [l]  In  experiment  the 
double  quantum  dot  is  formed  by  applying  appropriate 
voltages  to  external  gates  and  we  model  the  resulting 
confinement  potential  V2d(C)  by  a  square  well  poten¬ 
tial  that  is  —A/2  for  —a/2  —  L  <  (  <  —a/2,  A/2  for 
a/2  <  C  <  a/2  +  L,  and  Vg  otherwise.  The  length  of 
the  dots  is  L  and  a  is  the  width  of  the  interdot  barrier. 
As  discussed  previously  we  do  not  expect  our  results  to 
change  qualitatively,  if  a  smoother  potential  is  used.  In 
the  double-dot  system  at  finite  detuning,  the  depths  of 
the  dots  change  affecting  the  single  particle  energies.  In 
the  numerical  calculation,  we  determine  the  eigenspec- 
trum  of  Vi d  exactly.  The  main  effect  of  the  detuning 
is  an  energy  shift  ±A/2  to  the  single-particle  eigenener- 
gies,  where  “+”  corresponds  to  right  dot,  and  ”  to  the 
left  dot.  Using  Eq.  ©  and  neglecting  interdot  tunneling, 
the  energies  of  the  localized  left  and  right  single-particle 
orbitals  are  approximately 

ES  «  En(Vg  ±  A/2)  +  Ea  «  Ecn(Vg)  ±  A/2  +  Ea.  (7) 

When  more  than  one  electron  is  confined  in  the  single 
or  double  dot,  electron-electron  interactions  become  im¬ 
portant.  The  electrons  interact  through  the  long-range 
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Coulomb  potential 


Vc(r1,r2) 


kd |ri  -  r2|  ’ 


(8) 


where  kd  denotes  the  dielectric  constant.  Coulomb  inter¬ 
actions  allow  for  certain  off-diagonal  matrix  elements  in 
valley  space  that  are  produced  by  intervalley  scattering^!]. 
However,  these  matrix  elements  are  small  for  quantum 
dots  with  a  size  much  larger  than  the  interatomic  dis¬ 
tance;  they  are  neglected  in  this  work18. 

To  obtain  an  accurate  description  of  interacting  few 
electron  systems,  we  extend  the  single-dot  treatment  of 
Refs.  1 181 H)]  to  the  double  dot  system.  We  construct  sin¬ 
gle  particle  orbital  localized  in  the  left  and  right  dots 
from  the  exact  single  particle  solutions  of  the  double  dot 
and  we  used  them  to  expand  the  many-body  Hamiltonian 
(see  more  details  in  Appendix  [A|.  The  single  particle  or¬ 
bitals  have  a  weak  dependence  on  the  spin- valley  degrees 
of  freedom  that  comes  from  the  dependence  of  the  wave 
vectors  k  and  k  on  r  and  s.  This  leads  to  a  dependence 
of  the  interaction  matrix  elements  on  the  spin-valley  de¬ 
grees  of  freedom.  However,  this  dependence  is  very  weak 
and,  to  a  very  good  approximation,  can  be  neglected. 
Thus,  interactions  can  be  considered  spin- valley  indepen¬ 
dent  which  allows  a  separation  of  orbital  and  spin-valley 
contributions  of  the  two-electron  solution^1®. 


B.  Separating  orbital  from  spin- valley  degrees  of 
freedom 

Since  interactions  can  be  considered  diagonal  in  spin- 
valley  degrees  of  freedom,  orbital  and  spin-valley  part 
of  the  two-electron  solutions  provide  independent  contri¬ 
butions  to  the  energies  and  the  wave  functionf^.  The 
total  two-particle  wave  function  must  be  antisymmetric 
with  respect  to  particle  exchange  and  the  symmetry  of 
the  orbital  part  must  always  be  opposite  to  that  of  the 
spin- valley  part.  Thus,  the  two-particle  spectrum  can 
be  grouped  according  to  their  orbital  symmetry  or  par¬ 
ity  under  particle  exchange  in  multiplets  of  six  states  if 
the  orbital  part  is  symmetric  or  ten  states  if  the  orbital 
part  is  antisymmetric.  The  energy  splitting  between  dif¬ 
ferent  multiplets,  called  e,  depends  on  the  orbital  part, 
which  is  determined  by  electron-electron  interactions  and 
longitudinal  confinement,  and  it  is  generally  given  by  a 
correlated  state  that  is  represented  as  a  superposition  of 
various  two-electron  orbital  wave  functions.  The  energy 
relations  within  a  multiplet  are  exclusively  determined 
by  Ea  that  includes  the  orbital  and  spin  Zeeman  terms 
as  well  as  the  spin-orbit  coupling  Eq.  ©•  Therefore,  the 
spin-valley  part  of  the  wave  function  always  has  the  sim¬ 
ple  form  shown  in  Table  [3] 

The  magnetic  field  dependence  of  the  two  multiplets  is 
illustrated  in  Fig.  [3]  The  competition  between  spin-orbit 
coupling  and  orbital  Zeeman  energy  leads  to  a  ground 
state  crossings  in  the  multiplet  with  antisymmetric  spin- 
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FIG.  3:  Table:  spin-valley  multiplets,  here  |tisi,  t2s2)±  = 
(|tisi,t2s2>  ±  |t2s2, tisi))/\/2.  The  states  in  each  multiplet  are 
grouped  in  three  columns  according  to  their  spin-orbit  energy. 
Figures:  Schematic  magnetic-field  dependence  of  antisymmet¬ 
ric  (a)  and  symmetric  (b)  multiplets. 


valley  part  at  a  critical  magnetic  field  Bc  =  Aso/^Horb- 
For  B  <  Bc  the  ground  state  of  the  antisymmetric  spin- 
valley  multiplet  [Fig.  [3]  (a)]  is  in  a  superposition  of  spin 
singlet  and  triplet,  since  \K  \.,K'  |)_  =  (| KK')+  ||t)_  + 
| KK')~  |4_t)  +  )/v/2-  For  B  <  Bc ,  the  lowest  state  of  the 
antisymmetric  spin-valley  multiplet  |  K  f  ,K  J,)  is  anti¬ 
ferromagnetic  in  spin  but  ferromagnetic  (or  polarized) 
in  valley  space.  The  lowest  state  of  the  symmetric  spin- 
valley  multiplet  | K  |,  K  j.)+  is  ferromagnetic  in  both  spin 
and  valley  space  for  all  positive  magnetic  fields. 

The  energy  difference  between  the  different  multiplets  e 
and  the  magnetic  field  B  determines  the  spin-valley  sym¬ 
metry  of  the  ground  state.  The  thick  curves  in  Figs.  [3](a) 
and  (b)  correspond  to  the  possible  ground  states  which 
we  label  according  to  their  spin  and  valley  symmetry. 


C.  Model  Description 

One  of  the  objectives  of  this  study  is  to  be  able  to  de¬ 
scribe  the  evolution  of  the  spectrum  as  the  detuning  is 
changed  from  small  to  large  and  the  low  energy  configu¬ 
rations  change  from  \\L,  If?),  i.e.  one  electron  per  dot, 
to  \2L,0R)  or  \0L,2R),  two  electrons  in  the  same  dot. 
An  accurate  description  of  the  strong  interaction  effects 
requires  the  inclusion  of  several  single  particle  bands  of 
the  double  dot  system.  For  example,  the  behavior  of  the 
|2L,0i?)  configurations  is  expected  to  be  very  similar  to 
that  of  the  doubly  occupied  single  dot  and  in  that  sit¬ 
uation  the  strong  correlations  need  to  be  described  by 
many  single  particle  orbitals18  19 .  This  situation  makes 
the  description  and  the  interpretation  of  the  results  not 
very  intuitive.  However,  we  can  significantly  simplify  the 
description  by  realizing  that  no  matter  how  strongly  cor¬ 
related  the  system  is,  the  parity  is  a  good  quantum  num¬ 
ber  and  the  states  can  be  classified  according  to  the  state 
parity.  Thus,  we  can  model  the  exact  system  with  a  sim- 
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pie  effective  Hamiltonian  in  the  charge  degrees  of  freedom 
that  captures  this  dependence  on  parity  and  describes  the 
energetically  lowest  multiplets  of  the  \1L,  1R),  \2L,0R) 
and  |0L,  2i?)  configurations. 

The  charge  degrees  of  freedom  of  two-electrons  in  a 
double  dot  can  have  three  configurations:  \2L10R)±, 
\1L,  li?)  ,  and  \0L,  2i?)  ,  where  ±  characterizes  the 
conserved  orbital  symmetry  (i.e.  +  for  antisymmetric 

spin- valley  states  and  -  for  symmetric  spin- valley  states). 
Within  this  model,  interaction  effects  can  be  obtained  by 
diagonalizing  the  effective  Hamiltonians  Hg  and  Has , 
where  S  and  AS  denote  symmetric  and  antisymmetric 
spin-valley  state.  The  complex  multiband  problem  is 


then  reduced  to 

simple 

3  times  three  matrices. 
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These  effective  Hamiltonians  include  the  onsite  and  near¬ 
est  neighbor  interactions  V  and  Vlr,  the  tunnelings  in 
symmetric  ts  and  antisymmetric  tAS  configuration  and 
the  detuning  effects  A.  The  dependence  of  the  inter¬ 
action  on  the  symmetry  is  introduced  by  an  effective 
exchange  term  Vex  that  favors  the  antisymmetric  spin- 
valley  configuration.  Equations  [9]  and  [lO]  describe  the 
energy  related  with  the  orbital  part  of  the  wavefunction. 
The  total  energy  also  contains  the  contribution  of  the 
spin- valley  part  ESv  =  J2aiEi  (Vg)  +  Ea\na,  where  Ea 
was  defined  in  Eq.  ([6])  and  na  denotes  the  occupation  of 
states  with  spin-valley  a. 

To  gain  qualitative  understanding  of  the  interaction 
and  tunneling  terms  in  the  effective  Hamiltonian,  we  an¬ 
alyze  the  limiting  behaviors  of  the  low-energy  spectrum. 
First,  we  consider  the  limit  of  zero  detuning  and  large 
local  Coulomb  interactions  and  we  obtain  that  the  low¬ 
est  symmetric  and  antisymmetric  spin-valley  multiplets 
have  energies: 

^  »  e3v  +  vlx-v_*}°_v^+...,  (11) 

2t2 

Es  ~  Esv  +  Vlr  -  — —  s  +•••■  (12) 

V  -  Vlr  +  Vex 

In  the  strong  interaction  regime,  our  numerical  calcula¬ 
tions  indicate  that,  to  a  good  approximation  ts  ~  tAS  ~ 
t\2 .  This  approximation  allows  to  obtain  a  simple  expres¬ 
sion  for  the  energy  splitting  between  different  multiplets, 
e  =  E9S  —  E9a  ss  4 t\2Vex/V2  in  the  limit  of  zero  detuning 
and  for  VexCV-  VLR. 

For  a  biased  double  dot  system,  the  single  particle  en¬ 
ergies  acquire  an  energy  shift  of  ±A/2  and  in  the  limit 
of  large  detuning  the  two  electrons  occupy  the  same  dot. 


In  this  limit,  energies  of  the  lowest  symmetric  and  anti¬ 
symmetric  spin- valley  multiplets  are 

Ebas  «  Esv  —  A  +  V  —  Vex  +  ...,  (13) 

E%  »  Esv-  A  +  V  +  Vex  +  ....  (14) 

Thus,  the  energy  splitting  e  =  Eg  —  EA  S3  2Vex  in  the 
(2,0)  configuration  is  mainly  controlled  by  the  exchange 
mechanism. 

This  effective  model  allows  a  simple  and  intuitive  un¬ 
derstanding  of  the  underlying  physical  behavior  of  the 
double  dot  system.  However,  to  extract  the  parameters 
V,  Vlr ,  Vex  and  ti2  we  need  to  solve  exactly  the  single 
and  double  dot  system.  In  the  next  Section,  we  analyze 
the  behavior  of  the  double-dot  system  at  different  mag¬ 
netic  field,  detuning  and  interaction  regimes  by  compar¬ 
ing  the  exact  diagonalization  solutions  discussed  in  the 
previous  subsection  with  the  model  Hamiltonian.  From 
this  comparison,  we  extract  the  parameters  of  the  model. 


III.  RESULTS 

In  this  section,  we  analyze  the  energy  spectrum  of  two 
interacting  electrons  in  a  double  dot  in  three  different  de¬ 
tuning  regimes:  (A)  small  detuning,  with  one  electron  in 
each  dot;  (B)  strong  detuning,  with  both  electrons  in  the 
same  dot  and  (C)  at  the  crossover  between  both  regimes. 
In  subsection  |IIID[  we  study  the  transport  properties  of 
the  double  dot. 

The  analysis  of  the  spectrum  and  transport  properties 
is  characterized  the  energy  spacing  e  which  can  be  de¬ 
fined  in  the  exact  calculations  as  the  energy  difference 
between  the  lowest  energy  state  with  an  antisymmetric 
orbital  part  and  and  the  lowest  state  with  a  symmetric 
orbital  part  at  zero  magnetic  field  B  =  0.  In  all  detun¬ 
ing  regimes,  we  find  e  >  0,  in  agreement  with  the  Lieb 
Mattis  theorem27}  However,  Coulomb  correlations  can 
significantly  reduce  e. 

In  our  analysis,  we  have  considered  different  double 
dot  configurations  by  changing  the  length  and  depth 
of  the  dots  as  well  as  the  interdot  distance  and  found 
the  same  scaling  of  interaction  effects  as  discussed  in 
ReflTSl  In  this  section  we  present  results  for  the  pa¬ 
rameters  L  =  70  nm,  a  =  20  nm,  Vg  =  78  meV, 
R  =  2.5  nm.  These  parameter  lead  to  an  energy  splitting 
between  the  two  lowest  single-particle  longitudinal  modes 
of  hu>o  ~  10.6  meV  and  the  symmetric  double  dot  sup¬ 
ports  five  bound  states  per  quantum  dot  that  are  all  in¬ 
cluded  in  the  multi-band  calculation.  The  dielectric  con¬ 
stant  kd  is  varied  between  (1.5  <  kd  <  3.5)  which  allows 
us  to  explore  the  strongly  interacting  regime  where  new 
transitions  occurs.  Experimentally,  however,  it  is  easier 
to  change  the  length  of  the  dots.  The  parameter  that 
characterize  the  strength  of  the  interactions  is  the  ratio 
U/(hu>o)  between  typical  single  particle  excitation  ener¬ 
gies  Tuvq  and  the  typical  intradot  interaction  energy  scale 
U  =  e2/kdL.  This  ratio  is  typically  (1  <  U/(huJo)  <  5) 
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which  imply  a  moderate/strong  interaction  regime.  An¬ 
other  relevant  parameter  that  determines  the  regime  of 
a  quantum  dots  is,  as  we  will  discuss  below,  the  ratio 
2Vex/Aso  that  reflects  the  competition  of  strong  correla¬ 
tions  and  spin-orbit  effects  and  determines  the  spin- valley 
nature  of  the  ground  state. 


A.  Low  energy  spectrum  of  symmetric  double  dot 


A  [meV] 


Their  configurations  are  approximately 


1 M)  * 

*  IF,  1F)+  \K  4.;  K'  t)~  > 

(15) 

1  AfF)  * 

*  \IL,IR)+\K  ,  and 

(16) 

|FF)  * 

i  \1L,1R)-\K  UK  1} . 

(17) 

In  the  exact  diagonalization  calculation,  we  extract  the 
orbital  structure  and  conclude  that  |1  F,  1  R}±  is  given  to 
a  very  good  approximation  by  single  Slater  determinant 
formed  with  left  and  right  orbitals  in  the  lowest  band. 


FIG.  5:  (Color  online)  Parameters  of  the  effective  Hamilto¬ 
nian  as  a  function  of  the  dielectric  constant. 


FIG.  4:  (Color  online)  Low-energy  spectrum  of  a  double  dot 
as  a  function  of  B.  Solid  symbols  (red  online)  correspond  to 
states  belonging  to  the  antisymmetric  spin-valley  multiplet 
and  open  symbols  (blue  online)  to  the  symmetric  spin-valley 
multiplet.  The  solid  curves  correspond  to  the  effective  Hamil¬ 
tonian  description. 


First,  we  analyze  the  numerical  results  obtained  with 
the  multiband  treatment.  Figure [4]  shows  the  low-energy 
spectrum  as  a  function  of  B  for  kd  =  2.5  in  a  double  dot 
with  A  =  0  (zero  detuning) .  In  the  low-energy  spectrum, 
we  can  recognize  the  symmetric  and  antisymmetric  spin- 
valley  multiplets  discussed  in  Table/Fig.  [3]  The  energy 
difference  e  between  the  two  multiplets  is  not  observable 
in  the  energy  range  of  Fig.  [4]  Since  interdot  tunneling  is 
very  small  compared  with  the  interaction  energy,  the  two 
electrons  occupy  different  dots  to  avoid  strong  intradot 
interactions.  The  interdot  interaction  is  almost  indepen¬ 
dent  of  the  orbitals  occupied  in  each  dot.  Thus,  at  small 
detuning  there  is  a  negligible  occupation  of  higher  bands. 

In  Fig.  [4j  we  identify  three  states  that  will  be  rele¬ 
vant  for  the  discussion  of  Pauli  blockade.  One  of  them 
has  mixed  valley-spin  symmetry  and  we  label  it  as  |  M). 
The  second  state  is  antiferromagnetic  in  spin  and  ferro¬ 
magnetic  in  valley  degree  of  freedom  and  we  label  it  as 
| AfF).  The  third  state,  labeled  IFF),  is  ferromagnetic 
in  both  spin  and  valley  degrees  of  freedom.  States  |M), 
|  AfF)  belong  to  the  antisymmetric  spin- valley  multiplet 
and  |FF)  belongs  the  symmetric  spin- valley  multiplet. 


We  can  approximate  the  energies  of  |  M),  \  AfF)  and 
|FF)  using  Eqs.  (TFI),  (fTTb  and  (p~2]), 


«  2El{Vg)  ~2t\2/{V  -Vex)  +  VLR-  Asom 
E{aJf)  »  2Et{Vg)-2t\2/{V-Vex)  +  VLR-2Bnm 
Ew-f)  »  2Ec1(Vg)-2t212/(V  +  Vex)  +  VLR 

—A so  ~  2 B(fi  orb  T  H’spin)-  (20) 


At  B  =  0,  the  ground  state  | M)  is  only  separated  by 
the  very  small  superexchange  energy  e  ~  At\2Vex/V'2 
from  the  lowest  spin-valley  symmetric  triplet  that  con¬ 
tains  |FF).  These  four  states  are  separated  from  the  rest 
of  the  spectrum  by  a  much  larger  energy  scale  given  by 
A  so-  This  energy  structure  resembles  the  singlet-triplet 
splitting  in  GaAs  double  quantum  dots.  At  finite  fields, 
|FF)  is  the  ground  state.  The  first  excited  state  changes 
with  increasing  magnetic  field  from  |  M)  to  |A/F)  at 
Bc  =  ASo/(2Horb),  as  shown  in  Fig.  B  This  crossing  be¬ 
tween  two  antisymmetric  spin-valley  states  has  no  anal¬ 
ogous  in  standard  GaAs  quantum  dots. 

From  the  analysis  of  the  single  and  double  dot  spec¬ 
trum,  we  can  obtain  the  parameters  of  the  charge  ef¬ 
fective  Hamiltonian  (Eqs.  [9  10).  Figure  [5]  presents  the 
parameters  V,  Vrr,  H2  and  Vex  for  a  double  dot  with 


L  =  70  mn,  a  =  20  nm,  Vg  =  78  meV,  R  =  2.5  nm, 
and  A  =  0.  The  V  and  Vex  are  obtained  from  the  two- 
electron  spectrum  in  a  single  dot  and  Vlr  and  ti2  ob¬ 
tained  from  double  dot  spectrum.  The  black  solid  curves 
in  Fig.  [4]  show  the  prediction  from  the  model  using  the 
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parameters  from  Fig.  [5] 


B.  Low  energy  spectrum  for  large  detuning:  Two 
electrons  in  a  single  dot 


When  the  detuning  is  large  enough  to  overcome  the  in¬ 
tradot  interaction,  both  electrons  occupy  the  same  dot, 
and  the  charge  degree  of  freedom  of  the  low-energy  eigen¬ 
states  can  be  described  by  the  |2L,0i?)  configuration. 
Thus,  the  energy  spectrum  resembles  the  one  obtained 
for  two  electrons  in  an  isolated  dolP. 


0.0  0.1  0.2  0.3  0.4 

A  [meV] 


FIG.  6:  (Color  online)  Low  energy  spectrum  at  large  detun¬ 
ing  (A  ~  35  meV).  Solid  symbols  (red  online)  correspond  to 
states  belonging  to  the  antisymmetric  spin-valley  multiplet, 
and  open  symbols  (blue  online)  to  the  symmetric  spin-valley 
multiplet.  Curves  represent  the  model  predictions. 


Figure  [6]  presents  the  low  energy  spectra  for  large  de¬ 
tuning  for  kd  =  2.5.  The  solid  curves  represent  the  effec¬ 
tive  Hamiltonian  predictions.  The  multiband  structure 
of  the  solutions  introduce  corrections  to  the  spin-valley 
dependence  of  the  spectrum  which  can  be  absorbed  in  ef¬ 
fective  fxorb  and  Ago-  However,  this  corrections  are  small 
and  correspond  to  a  few  percent  changes  to  the  bare  fiorb 
and  Ago  values. 

At  zero  magnetic  field,  the  antisymmetric  spin-valley 
multiplet  is  favored.  This  is  in  agreement  with  Lieb  Mat- 
tis  theorem27  that  states  that  the  two  particle  ground 
state  always  has  a  symmetric  orbital  part.  In  our  two- 
electron  system,  we  can  understand  this  prediction  from 
the  analysis  of  the  orbital  symmetry  of  the  wave  function. 
In  the  noninteracting  limit,  the  orbital  ground  state  is 
constructed  with  both  electrons  in  the  lowest  band  cor¬ 
responding  to  a  symmetric  orbital  wave  function,  i.e. ,  an 
antisymmetric  spin-valley  wave  function.  To  form  an  an¬ 
tisymmetric  orbital  wave  function  at  least  one  electron 
has  to  occupy  an  excited  state.  Thus,  in  the  noninter¬ 
acting  limit,  e  is  given  by  the  level  splitting  hujo  between 
the  first  two  bands.  However,  the  energy  difference  e  be¬ 
tween  the  two  multiplets  shown  in  Fig.  [6] is  much  smaller 
than  huiQ.  Furthermore,  the  reduction  of  e  with  respect 
to  the  level  splitting  is  much  stronger  for  smaller  dielec¬ 


tric  constant  or  longer  dots,  since  interaction  effects  are 
more  relevant.  In  the  limit  of  infinite  interactions,  the 
electrons  are  strongly  localized  at  the  positions  that  min¬ 
imize  the  interaction  energy.  Then  the  orbital  symmetry 
of  the  wave  function  becomes  irrelevant  and  e  vanishes. 
This  effect  is  a  signature  of  a  formation  of  a  Wigner 
molecul^^SSl. 

In  the  effective  Hamiltonian,  the  formation  of  a  Wigner 
molecule  is  manifested  in  the  reduction  of  the  parameter 
Vex.  This  is  evident  in  Fig.  [5]  that  shows  that  Vex  is  a 
growing  function  of  the  dielectric  constant.  Even  though 
the  results  of  Fig.[5]are  for  A  =  0,  we  note  that  V  and  Vex 
depend  weakly  on  the  detuning  and  can  be  approximated 
by  A  =  0  predictions  for  all  detunings  studied  here  A  < 
35  meV.  In  contrast,  the  tunneling  t\2  is  strongly  affected 
by  the  detuning  and  is  reduced  almost  by  approximately 
a  factor  of  two  in  comparison  with  the  A  =  0  case. 

The  reduction  of  e  implies  that  the  exact  eigenstates 
become  strongly  correlated  and  cannot  be  written  as  non¬ 
interacting  wave  functions  to  characterize  them.  How¬ 
ever,  we  can  still  label  the  low-lying  (2,0)  states  accord¬ 
ing  to  their  conserved  quantum  numbers: 


1 M)  = 

\2L,0R)+\KUK'  tr, 

(21) 

1  AfF)  = 

\2L,  0R)  +  \K  J,;  I\  t)-  ,  and 

(22) 

|FF)  = 

\2L,  0R)~  \K  UK  1} 

(23) 

Here  the  states  |2Zv,  0-R)=t  are  correlated  orbital  states  of 
two  electrons  in  the  left  dot.  Our  numerical  calculations 
shows  that  several  bands  are  needed  to  represent  these 
states  accurately.  We  note  that  Secchi  and  Rontani  found 
qualitatively  the  same  correlation  effects  for  a  parabolic 
well  with  weak  confinement®,  which  suggests  the  robust¬ 
ness  of  these  correlation  effects.  We  note  that  because  of 
the  conserved  symmetries,  the  eigenstates  at  small  de¬ 
tuning  | M),  | AfF),  and  \FF)  will  evolve  in  | M),  \ AfF), 
and  |F.F)  for  large  detuning. 


Using  the  effective  Hamiltonian  along  with  the  approx¬ 
imate  description  of  the  single  particle  energies  [Eqs(7j 


13  and  14  ,  we  can  obtain  simple  expression  for  the  (2,0) 


configuration  energies: 


E\m)  ~  2FJ(US)  +  V  -  Uex  -  Aso  -  A,  (24) 

E\AfF)  ~  2 Ei(Vg)  +  V  —  Vex  —  2 Bfj,orb  —  A,  an<(125) 

E\ff)  ~  2Fj(V),)  +  V  +  Vex  —  Aso 

2H(/rorf)  A.  (26) 

At  B  =  0,  | M)  is  the  ground  state.  Above  a  critical 
magnetic  field,  there  is  a  ground-state  transition  to  ei¬ 
ther  |  AfF)  (  Ag0  <  e)  or  |FF)  (  A  so  >  e)  due  to  the 
orbital  Zeeman  term.  Thus,  the  reduction  of  e  leads  to  a 
ground  state  transition  to  the  ferromagnetic  state  |FF) 
at  finite  magnetic  field.  In  practice,  the  formation  of  a 
ferromagnetic  ground  state  can  be  experimentally  con¬ 
trolled  by  changing  the  length  of  the  dots. 


C.  Transition  from  a  double-  to  a  single-dot 
regime. 

We  now  analyze  the  transition  between  the  limiting 
behaviors  discussed  in  the  previous  two  sections. 


FIG.  8:  (Color  online)  e  as  a  function  of  detuning,  symbols 
correspond  to  the  numerical  results  and  the  curve  is  the  ef¬ 
fective  Hamiltonian  prediction. 


diction  as  shown  in  Fig.  [8j 


FIG.  7:  (Color  online)  Zoom  of  first  crossing  of  (1,1)  and  (2,0) 
states  for  B  =  0.  Solid  symbols  (red  online)  correspond  to 
states  belonging  to  the  antisymmetric  spin-valley  multiplet, 
and  open  symbols  (blue  online)  to  the  symmetric  spin-valley 
multiplet.  Curves  represent  the  model  predictions. 

Figure  [7]  shows  the  crossover  between  |lL,li?)  states 
and  |2L,  0i?)  states  of  the  double  dot  system  at  B  =  0. 
In  absence  of  interdot  tunneling,  the  two-electron  spec¬ 
trum  shows  sharp  crossings  between  the  |1L,  li?)  states 
and  |2L,0f?)  states  with  increasing  detuning.  Because 
of  the  interdot  tunneling,  crossings  between  states  with 
the  same  symmetries  turn  into  avoided  crossings.  The 
avoided  crossings  occur  at  the  same  critical  detuning 
for  states  of  the  same  multiplet.  The  avoided  crossings 
within  the  multiplet  with  an  antisymmetric  orbital  part 
occur  at  relatively  larger  detuning  since  the  tunneling 
electron  is  forced  to  occupy  an  excited  band. 

The  crossover  regime,  presented  in  Fig.  [7J  is  strongly 
affected  by  interactions.  The  difference  between  the  criti¬ 
cal  detunings  belonging  to  the  avoided  crossings  of  states 
with  symmetric  and  antisymmetric  orbital  part  is  a  di¬ 
rect  measure  for  the  energy  splitting  e.  Furthermore, 
correlations  in  the  \2L,  OR)  states  decrease  the  tunneling 
coupling  to  the  corresponding  |1  L,  If?)  state,  leading  to 
a  sharper  avoided  crossing. 

In  the  effective  Hamiltonian  description,  the  lowest 
\IL,IR)  and  \2L,  OR)  configurations  are  close  in  energy 
when  A  ~  V  —  Vlr-  In  this  regime,  the  |0L,  2 R)  con¬ 
figurations  are  energetically  suppressed  and  do  not  affect 
the  low  energy  spectrum.  Within  this  approximation,  we 
can  obtain  a  simple  expression  for  e, 

e  =  Vex  +  2  4tj2  +  (A  —  Vex  —  V  +  Vlr)2 

-  V^4 i?2  +  (A  +  Vex  -  v  +  VLR)^J  (27) 
This  expression  compares  well  with  the  numerical  pre- 


FIG.  9:  (Color  online)  Phase  diagram  for  a  double  quantum 
dot  with  five  bands  and  kd  =  1.5  (a)  and  kd  =  3.5  (b).  Black 
solid  lines  separate  regions  with  different  spin- valley  symme¬ 
try.  Regions  with  the  same  greyscale  (  color  online)  charac¬ 
terize  the  ground  state. 


From  the  analysis  of  the  low-energy  spectrum,  e.  g., 
Fig.0  we  can  extract  the  phase  diagram  as  a  function 
of  the  detuning  A  and  the  magnetic  field  B.  The  phase 
diagrams  for  strong  (kd  =  1.5)  and  weak  (kd  =  3.5)  inter¬ 
actions,  presented  in  Fig.  |9j  exhibit  clear  differences.  At 
exactly  zero  magnetic  field,  the  ground  state  always  has 
a  symmetric  orbital  wave  function  (| M)  and  | M)).  This 
state  is  preferred  with  respect  to  other  states  within  the 
same  multiplet  by  the  spin-orbit  coupling  A  so  and,  with 
respect  to  states  of  the  other  multiplet,  by  the  energy  gap 
e.  The  ground  state  changes  at  a  critical  magnetic  field 
to  a  valley-polarized  state.  The  valley  polarized  ground 
state  has  symmetric  orbital  part  for  A  so  <  £  and  an  an¬ 
tisymmetric  orbital  part  for  A  so  >  £■  For  small  detuning 
A,  the  energy  splitting  e  is  caused  by  the  small  superex¬ 
change  and  e  <  A  so-  However,  for  large  detuning,  both 
electrons  are  on  the  same  dot,  and  the  value  of  e  can 
be  larger  or  smaller  than  A  so-  Thus,  if  2Vex  >  A  so, 
there  is  a  transition  to  a  valley-polarized  state  with  a 
symmetric  orbital  part,  | AfF),  as  detuning  is  increased 
[as  observed  in  Fig.  IU(b)]. 


9 


D.  Sequential  transport 


In  this  section,  we  discuss  how  the  two-electron  eigen- 
spectrum  affects  sequential  transport  through  a  double 
dot.  In  particular,  we  analyze  the  existence  of  Pauli 
blockade  in  serially  coupled  dots  that  would  lead  to  a 
current  rectification  in  DC  transport6  Figure [lC)| shows 
a  schematic  description  of  two  situations  when  a  Pauli 
blockade  occurs.  In  both  cases,  the  double-dot  system  is 
detuned  such  that  there  is  always  at  least  one  electron 
in  the  left  dot.  States  relevant  for  transport  through  the 
double  dot  are  the  |2L,0I?)  states  | M)  and  | FF)  as  well 
as  the  \1L,  1 R)  states  | M)  and  | FF).  In  Fig. 
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the  ver¬ 
tical  axes  denotes  energy  and  the  horizontal  axes  denotes 
the  spatial  coordinate  along  the  nanotube.  The  center  of 
each  figure  shows  the  double  dot  created  by  tunnel  bar¬ 
riers  to  the  contacts  and  between  the  dots.  The  gray 
rectangles  to  the  left  and  right  of  the  double  dot  are  the 
Fermi  seas  in  the  contacts. 


(a)  B=0  Blocked  (b)  B=B*  Blocked 


FIG.  10:  (Color  online)  Schematic  representation  of  a  Pauli 
blockade  for  zero  magnetic  field  and  negative  bias  (a),  and 
finite  magnetic  field  and  positive  bias  (b).  The  white  areas 
in  the  center  of  each  figure  represent  the  two  dots.  The  black 
circle  in  the  left  dot  represent  the  position  and  energy  of  one 
of  the  electrons.  The  states  |  M)  and  |  M)  are  represented  by 
solid  lines  with  | FF)  and  | FF)  by  dashed  lines.  These  states 
are  organized  vertically  according  to  their  energies  and  also 
represent  the  possible  positions  of  the  second  electron.  The 
wide  gray  rectangles  next  to  the  left  and  right  dots  represent 
the  energy  bias. 


For  a  positive  bias  voltage  (opposite  to  Fig.  10  (a))  ,  the 
electrochemical  potential  in  the  left  contact  is  larger  than 
in  the  right  one,  and  current  flows  via  a  sequential  tunnel¬ 
ing  process  |1L,0I?)  — >  \2L,  OR)  —$■  \1L,IR)  — >  |1L,0I?). 
Interdot  tunneling  is  assumed  to  conserve  the  spin- valley 
degree  of  freedom  and  allows  for  transitions  between 
states  | M)  and  | M)  or  | FF)  and  | FF).  Because  of  the  se¬ 
quential  transport  setup,  the  left  dot  couples  to  the  left 
reservoir,  allowing  for  transitions  between  the  \1L,  OR) 
and  |2L,0f?)  states.  Analogously,  the  right  contact  al¬ 
lows  for  transitions  between  the  |1L,  01?)  and  \1L,1R) 
states.  Figure  10  (a)  considers  Pauli  blockade  at  B  =  0. 

In  this  scenario,  the  current  is  blocked  when  state  | FF) 
is  occupied  by  an  electron  tunneling  in  from  the  right 
reservoir.  In  fact  a  B  =  0  the  ground  state  of  the  sym¬ 
metric  multiplet  is  three  fold  degenerate,  but  as  long 


as  spin  and  valley  are  conserved  this  degeneracy  does 
not  change  the  physics,  see  Appendix  [B]  Once  in  state 
| FF),  the  electron  cannot  tunnel  back  because  of  the 
filled  Fermi  sea  in  the  right  contact.  It  also  cannot  tun¬ 
nel  to  the  left  dot  since  the  corresponding  state  |  FF) 
has  much  larger  energy  than  state  |  FF)  because  of  the 
energy  difference  e  between  states  | M)  and  IFF).  How¬ 
ever,  e  can  be  strongly  reduced,  allowing  for  a  finite  exit 


rate  from  state  |FF)  as  shown  in  Fig.  10  (a).  This  effect 


might  explain  the  absence  of  a  Pauli  blockade  in  Ref.  [TT] 
This  Pauli  blockade  implies  a  rectification  of  the  current 
since,  by  inverting  the  bias  voltage,  a  finite  current  can 
flow. 


In  order  to  make  these  statements  more  quantitatively 
we  calculate  the  stationary  current  with  a  rate  equa¬ 
tion  approaclP^So].  We  describe  the  regime  of  possible 
blockade  depicted  in  Fig[l0](a)  while  assuming  spin-orbit 
coupling  to  be  much  bigger  than  temperature,  transport 
voltage  and  external  coupling  the  the  reservoirs. 

The  rate  equations  and  the  resulting  stationary  current 
are  given  in  Appendix  [B]  The  dependence  of  the  station¬ 
ary  current  Ion  in  the  blockade  setup  and  the  probability 
to  be  in  any  of  the  three  degenerate  |FF)  states  are  given 
by: 


J 

B  +  C{2t%  +  e2) 

(28) 

C{tl  +  e>)  +  D 

FF  B  +  C{2t2s  +  e2) 

(29) 

Here  e  denotes  the  energy  splitting  between  the  states 
|  M)  and  |FF)  and  ts  is  the  interdot  tunneling  rate  be¬ 
tween  states  |FF)  and  |FF).  The  constants  A,B,C,D 
depend  on  the  coupling  strength  to  the  contacts  and  spec¬ 
tral  weights  for  the  tunneling  probabilities  and  their  form 
is  given  in  Appendix  |B| 

If  e,  is  the  dominant  energy  scale  then  the  current  is 
suppressed  as  Ion  oc  1/e2  and  the  double  dot  gets  stuck 
in  states  |FF)  with  PfF  — >  1.  This  is  the  regime  of 
Pauli-blockade.  However,  we  find  that  due  to  interaction 
effects  interdot  tunneling  ts  can  even  exceed  the  energy 
splitting  e,  thus  removing  the  blockade  mechanism. 


Due  to  Pauli-blockade  the  double  dot  acts  as  a  cur¬ 
rent  rectifier,  since  by  reversing  the  transport  voltage  in 
Fig.  [To]  a)  a  finite  current  can  flow.  Fig.  11  illustrates 
how  effective  the  double  dot  acts  as  current  rectifier  as 
function  of  e. 


At  finite  magnetic  field,  the  transport  behavior  can 
also  be  very  interesting.  For  example,  if  interactions  are 
strong  enough  to  suppress  e  to  become  on  the  order  of 
the  spin-orbit  coupling  (A so  >  e)  but  still  much  larger 
than  the  interdot  tunneling  (e>f),  then  a  current  block¬ 
ade  can  occur  in  the  opposite  bias  direction  of  the  zero 
field  case.  This  situation  is  depicted  in  Fig.  10  (b).  A 
magnetic  field  B*  ss  e/2(p,orb  +  iispin)  is  applied  so  that 
states  |FF)  and  | M)  are  degenerate,  while  in  the  |1L,  li?) 
charge  configuration,  the  state  |FF)  is  ground  state  and 
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kd 

1.5  2.  2.5  3.  3.5 


FIG.  11:  (Color  online)  Ratio  between  leakage  current  (28 1 
corresponding  to  setup  in  FiaTola)  and  the  finite  current  for 
reversed  transport  voltage  (|B3[).  Used  parameters  Fl  =  F_r  = 
O.OlmeV,  tAS  =  ts  =  O.lmeV  and  Sm  =  0.7  Sf  =  1. 


separated  from  state  |  M)  by  the  energy  e  (because  of  the 
applied  field  B*).  The  physical  situation  of  Figure [lO] (b) 
resembles  that  of  Fig.  [lO]  (a),  and  also  results  in  a  cur¬ 
rent  blockade.  The  ability  to  control  the  direction  of  the 
current  rectification  by  varying  the  magnetic  field  can 
have  applications  to  carbon  nanotube-based  spintronic 
proposals. 


E.  ID-Quantum  dot  arrays 

The  rich  physics  associated  with  the  strong  correla¬ 
tions  of  double  occupied  dots  is  expected  to  have  sig¬ 
nificant  impact  in  the  many-body  behavior  of  a  one¬ 
dimensional  chain  of  dots.  The  double  dot  analysis  car¬ 
ried  out  here  can  be  used  as  a  starting  point  to  study  the 
behavior  of  a  linear  array  of  coupled  dots.  According  to 
our  analysis,  there  exists  an  interaction  regime  in  which 
the  behavior  of  two  electrons  in  the  same  dot  is  strongly 
correlated  and  controlled  by  e  while  the  interaction  be¬ 
tween  electrons  in  different  dots  does  not  introduce  band 
mixing.  In  order  to  generate  an  effective  Hamiltonian 
that  captures  the  essential  effects  of  strong  onsite  corre¬ 
lations,  we  apply  the  ideas  of  Hubbard  operators  to  de¬ 
scribe  double  occupied  dots  or  doublons  (see  i.e.  Ref.  EH)- 
This  description  assumes  that  the  strong  onsite  Coulomb 
repulsion  suppresses  the  probability  to  occupy  a  dot  with 
more  than  two  electrons.  This  is  a  good  approximation 
for  low  enough  fillings. 

Formally,  we  start  from  the  complete  Hamiltonian  that 
describes  a  chain  of  dots.  This  Hamiltonian  represents 
an  extension  of  the  double  dot  Hamiltonian  introduced 
in  the  Appendix. 

We  introduce  the  operators  a£  n  a  which  create  an  elec¬ 
tron  in  the  r  dot,  in  the  n  orbital  state  and  with  a  given 
spin-valley  configuration  a.  To  simplify  the  notation  we 
introduce  a£  a  =  a'r  i  a  for  electrons  in  the  lowest  band 
which  describe  well  the  singly  occupied  dots.  The  state 


that  describes  two  electrons  in  the  r  dot  can  be  expanded 
in  single  particle  basis  as 


\dri,r,a,a')  —  fa,a'  ^ 


,t 


77,71, m  \u'r,n,au'r,m,a' 


n<.m 


+  (-1)X,n,a'ar,m,a)l°) 


(30) 


Here,  77  labels  the  spin-valley  symmetry  of  the  two- 
electron  state  77  =  1  for  antisymmetric  spin-valley  states 
and  77  =  2  for  symmetric  spin-valley  states  and  fa,a'  = 
yjl  +  Sa^a'.  Equation  30  defines  the  doublon  opera¬ 
tor  such  that  I  =  4raa'l°)  &Ild  dl  r  a  a'  = 


(— l)vd'  r  ai  a-  Note  that  the  Pn,n,m  do  not  depend  on 
a  and  a'.  The  coefficients  /3v,n,m  can  be  obtained  by 
diagonalizing  the  local  part  of  the  Hamiltonian  which 
amounts  to  solve  two  electrons  in  a  single  dot.  We  now 
assume  that  the  occupation  of  higher  excited  two-particle 
states,  which  are  separated  by  an  energy  gap  of  the  or¬ 
der  of  the  single  particle  level  spacing,  can  be  neglected. 
This  should  be  a  good  approximation  for  the  low  energy 
spectrum. 


The  effective  Hamiltonian  is  obtained  by  projecting  the 
complete  Hamiltonian  onto  the  subspace  of  empty,  singly 
and  doubly  occupied  dots  and  can  be  written  as 


Hef  f  =  P(He  +Hd  +  Hed  +  VC)P.  (31) 


where  P  =  Pi  represent  a  projector  to  the  physi¬ 
cally  allowed  subspace  where  Pi  =  |0)i(0|  +  ]CQ  \a)i(a\  + 
J2n  \dv)i{dr1\  (see  Ref.  [32]).  He  describes  the  behavior  of 
single  occupied  dots,  Hd  describes  the  behavior  of  dou¬ 
ble  occupied,  Hed  contains  a  coupling  between  singly  and 
doubly  occupied  dots  and  Vc  represents  the  Coulomb  in¬ 
teraction  between  different  sites.  The  explicit  form  of 
these  contribution  is: 


He  = 
Hd  = 

Hed  = 


^  '  Af-.r*  c)Yy  Urey  d-  ^  '  ir^'^roflr'a 
r,Oi  (r,r'),a 

E2b,1]  d)  ,d 

/  j  r,a,a'  r},raa'  r],raa' 

r,r),oi<ai' 

1  V  t{d)  A 

'  /  v  °r},r,r'Ujr],raa,u'r],r' aa' 

(r,r') 

^  '  ^r,r'  faia3fa2a3dllraia3ar,a2araidri 


(32) 


(33) 


(r,r)  ,r),r)  , 

Oil,  CK2,CK3 


^  ^  f aot' 9r),r,r'  \^Jr},r,OiOi,^Jr->°i^Jr' i0-'  H*c] 

(r,r')r],ct,cx' 


Vc  =  UeY, 


r^r' 


nrnrf 
r  —  r' 


(34) 


The  diagonal  part  of  He  and  Hd  represent  the  elec¬ 
tron  and  doublons  energies  composed  of  the  spin-valley 
energy  and  onsite  energies  (Er,a  =  Ef(Vr)  +  Ea  and 
E2bava,  =  2E^(yr)+Ea  +  Ea,+V  +  (-l)Wex)  ,  and  t  and 
the  tunnelings  of  the  electrons  and  doublons.  The 
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Hed  contains  a  term  that  describes  the  hopping  of  an  elec¬ 
tron  from  the  doubly  occupied  site  to  a  singly  occupied 
site  which  in  this  language  corresponds  to  an  exchange 
of  a  doublon  and  an  electron.  The  second  term  in  Hed 
represents  the  hopping  of  an  electron  from  a  singly  occu¬ 
pied  site  to  a  singly  occupied  site  and  viceversa,  which  is 
represented  as  destroying  two  atoms  and  creating  a  dou¬ 
blon.  Finally,  we  have  the  offsite  Coulomb  interaction  Vc 
in  terms  of  the  dot  density  nr  defined  as 

fir  ^  '  Ura(lra  +  2  'y  '  ^raa'^TOo:''  (^) 

Here,  the  summation  {a,  a'}  is  restricted  to  a  <  a!  for 
??  =  1  and  a  <  a1  for  ??  =  2.  For  the  offsite  interaction 
we  assume  exclusively  capacitive  coupling  and  neglect  a 
dependence  on  the  symmetry  of  the  two  particle  states. 

In  the  most  general  case  the  parameters  that  describe 
the  effective  Hamiltonian  depend  on  lattice  sites  posi¬ 
tions  and  can  be  controlled  by  changing  the  detuning  in 
each  lattice  site  Ar.  For  example,  an  enhancement  of  the 
superexchange  interactions  can  be  achieved  by  detuning 
some  of  the  dots  and,  therefore,  reducing  the  energy  cost 
of  double  occupancies. 

The  effective  Hamiltonian  can  be  applied  to  describe 
an  array  of  coupled  dots  in  many  different  regimes  and 
its  phase  diagram  might  exhibit  novel  phases.  In  par¬ 
ticular,  the  existence  of  spin-valley  “triplet”  close  in  en¬ 
ergy  to  the  spin-valley  “singlet”  can  lead  to  phenomena 
richer  in  comparison  to  the  standard  single-band  Hub¬ 
bard  model33.  For  example,  for  one  particle  per  site  (fill¬ 
ing  n  =  1),  it  is  known  that  the  ground  state  of  the 
usual  Hubbard  model  has  infinite  susceptibility  to  spin 
dimerization^,  i.e.  formation  of  singlet /triplet  bonds  be¬ 
tween  nearest-neighbor  sites.  However,  the  ground  state 
is  not  dimerized  since  the  formation  of  singlet  bonds,  say 
at  sites  (2  i,  2  i+1),  is  penalized  by  the  large  energy  cost 
of  the  remainder  triplet  components  between  sites  (2  i+1, 
2  i+2).  In  carbon  nanotubes  operating  in  the  strongly  in¬ 
teracting  regime  ,  e<  1,  even  spin- valley  “triplet”  states 
can  lower  their  energy  by  virtual  hoping.  The  later  sit¬ 
uation  reduces  the  energy  cost  of  triplet  formation  and 
the  infinite  susceptibilities  to  spin  dimerization  might  in 
this  case  translate  in  actual  dimerization  of  the  ground 
state. 

Away  from  n  =  1,  the  presence  of  a  spin- valley 
“triplet”  can  also  significantly  affect  the  magnetic  struc¬ 
ture  of  the  system.  In  particular,  for  systems  which 
already  exhibit  ferromagnetism  in  the  standard  single¬ 
band  Hubbard  model,  the  inclusion  of  the  spin-valley 
“triplet”  can  strengthen  and  extend  the  ferromagnetic 
phase.  For  example,  generic  2  and  3  dimensional  square 
lattice  geometries,  and  others  with  similar  connectivity 
conditions,  exhibit  Nagaoka  ferromagnetisnpS  when  the 
gain  of  kinetic  energy  of  a  single  hole  exceeds  the  decrease 
of  superexchange  energy.  The  latter  condition  is  fulfilled 
at  very  large  on-site  interactions.  In  quantum  dots  in 
carbon  nanotubes  with  small  e,  Nagaoka-type  ferromag¬ 


netism  might  become  stable  at  reduced  value  of  the  inter¬ 
action  since  the  virtual  hopping  of  the  spin- valley  triplets 
reduces  the  super-exchange  penalty  of  having  a  polarized 
state. 

For  arrays  with  filling  1  <  n  <  2  and  zero  detuning,  the 
low  energy  physics  consists  of  doubly  and  singly  occupied 
dots  while  the  occurrence  of  empty  sites  is  strongly  sup¬ 
pressed  since  it  implies  an  increase  of  double  occupancies. 
This  implies  that  the  number  of  electrons  and  doublons 
will  be  independently  conserved  and  that  the  second  term 
of  Hed  can  be  neglected.  Also,  the  tunneling  of  the  elec¬ 
trons  and  the  doublons  can  be  neglected.  Thus,  the  only 
relevant  terms  in  the  effective  Hamiltonian  are  the  on¬ 
site  energies,  the  long-range  Coulomb  and  the  electron- 
doublon  exchange  (second  term  in  Hed).  This  regime  can 
lead  to  interesting  phenomena  when  the  electron-doublon 
exchange  becomes  comparable  the  multiplet  splitting  e. 

Finally,  it  should  be  pointed  out  that  the  long-range 
Coulomb  interaction  can  have  a  preponderant  influence 
in  the  charge  distribution  in  small  arrays  of  quantum 
dots. 


IV.  SUMMARY  AND  CONCLUSIONS 

We  have  presented  a  detailed  study  of  the  few-electron 
eigenspectrum  of  a  double  quantum  dot  in  a  semicon¬ 
ducting  carbon  nanotube.  We  showed  how  the  spin- 
valley  physics  leads  to  the  formation  of  multiplets.  The 
internal  energy  structure  of  the  multiplet  is  practically 
unaffected  by  a  change  in  either  the  confinement  poten¬ 
tial  or  the  interaction  strength,  but  the  energy  gap  be¬ 
tween  different  multiplets  is  strongly  modified  by  both. 
We  showed  that  for  sufficiently  strong  interactions,  the 
spin-orbit  coupling  can  exceed  the  energy  splitting  be¬ 
tween  states  with  symmetric  and  antisymmetric  orbital 
parts  for  any  detuning  between  the  dots.  This  situa¬ 
tion  modifies  the  two  particle  phase  diagram.  Above  a 
critical  interaction  strength,  the  ground  state  at  small,  fi¬ 
nite  magnetic  fields  is  always  ferromagnetic  independent 
of  the  detuning.  Furthermore,  in  this  strongly  interact¬ 
ing  regime,  the  blockade  of  linear  transport  gradually 
disappears  since  the  reduction  of  the  multiplets’  energy 
splitting  e  allows  a  finite  tunneling  probability  from  the 
|1L,  II?)  state  to  the  |2L,0I?)  state  even  for  a  symmet¬ 
ric  spin- valley  part.  Pauli-blockade  physics  will  occur  for 
weak  enough  Coulomb  correlations,  which  can  be  sup¬ 
pressed  by  either  working  with  short  dots,  or  by  covering 
the  nanotube  by  strong  dielectric^.  We  note  that  a  well 
developed  Pauli-blockade  is  the  precondition  for  realiza¬ 
tion  of  spin-qubits  in  double  dots  allowing  for  coherent 
singlet-triplet  manipulatiorPESl. 

Our  understanding  of  the  double  dot  physics  can  be 
used  as  an  starting  point  to  analyze  the  behavior  of  an 
array  of  coupled  dots.  The  effective  Hamiltonian,  as  pre¬ 
sented  in  Eq.  [3TJ  can  be  applied  to  describe  an  array  of 
coupled  dots  in  certain  regimes.  In  the  future,  we  would 
like  to  explore  the  many-body  physics  associated  with 
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the  strong  onsite  correlations  and  its  consequences  in  the 
magnetic  ordering  of  the  system. 

In  the  analysis  presented  here,  we  neglected  terms  that 
flip  either  the  spin  or  valley  degree  of  freedom,  or  both. 
However,  such  mechanisms  could  change  the  various  level 
crossings  to  avoided  crossings  and  thus  open  new  ways  to 
control  the  spin-valley  degree  of  freedom  as  well  as  the 
transport  properties  of  this  systerrP^U. 
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Appendix  A:  Coulomb  matrix  elements 


If  the  rest  of  the  integrand  does  not  depend  on  <pi  and 
ip2,  the  integration  over  those  coordinates  can  be  done 
analytically, 


^  /»27T  /*27T 

I  I  K0(R|2gsin[(<£i  —  ip2)/2]\)dip1dip2  = 

Io(RM)K0(R|<z|).  (A5) 

Also,  the  integration  over  </i  and  £2  can  be  done  ana¬ 
lytically  since  $(£)  has  a  simple  functional  dependence 
that  can  be  written  in  terms  of  exponentials  with  real 
or  imaginary  arguments.  Once  those  integrals  are  ob¬ 
tained,  we  are  left  with  the  momentum  integration  over 
q  in  Eq.  (A4).  Defining 


7nP(q)  =  I  eiq(l  {$t  (Ci)$p(Ci)}  dCi,  (A6) 


we  can  write  the  Coulomb  interaction  matrix  as  a  one- 
dimensional  numerical  integration  over  the  momentum 

q, 


We  create  an  orthonormal  set  of  localized  single  parti¬ 
cle  orbitals  <!>,  solving  exactly  the  double  dot  system  and 
we  use  this  set  to  construct  the  many-body  Hamiltonian 
for  the  interacting  double-dot  system: 


H  =  Y^Ei  a\ai+Y^  Tij  a\ar 
i 


\  y~]  Uijk{a\a^akai- 

ijkl 


(Al) 


Here  E,  =  are  the  single-particle  energies, 

Tjj  =  (|l>j|.£/o|<hj)  are  the  tunneling  matrix  elements, 
and  Uijki  =  correspond  to  the  interac¬ 

tion  matrix  elements.  At  zero  detuning,  the  tunneling 
matrix  elements,  T^,  are  only  nonzero  between  states 
with  the  same  <7,  r,  n  and  located  at  different  wells.  For 
nonzero  detuning,  tunneling  between  different  bands  be¬ 
come  possible.  The  electrons  interact  through  the  long- 
range  Coulomb  potential,  Eq.  [8j  and  the  interaction  ma¬ 
trix  elements  take  the  form 


Un 


47 T2k, 


J  <(Ci)  •  SP(Ci) 


C(C2)-<MC2)-#idCiWC2,  (A2) 


Un 


e 

kdTT 


Fnp(q)Fmr(-q)  l0{R\q\)K0{R\q\)dq. 

(A7) 


In  the  calculation  of  the  Coulomb  matrix  elements,  we 
neglect  the  dependence  of  the  longitudinal  wavevec- 
tor  kn  on  spin,  valley,  and  magnetic  field.  As  dis¬ 
cussed  in  Sec.  m  this  approximation  is  good  for  k 
®ab/($qR),  Ago/( 2tiv),  which  is  certainly  true  for  the 
parameters  used  in  this  paper. 


There  is  another  effective  interaction  that  takes  into 
account  intervalley  scattering: 


U«£r  =  ^j  *t(Cl)  '  *p(Cl) 

C(C2)  •  M(2)VKK' dvidCickpzdb  (A8) 

where  VKR  is  a  short  range  interaction  that  can  be  taken 
proportional  to  <5(£i  —  </2)-  We  have  estimate  the  URRpr 
terms  and  we  have  concluded  that  they  only  introduces 
minor  corrections.  For  that  reason,  we  neglect  this  con¬ 
tribution  in  the  present  work. 


Appendix  B:  Rate  equations 


where 

i  _  i  1 

r  !ri  ^*2 1  y/2R2[l  -  cos(^i  -  ip2)]  +  (Ci  -  C2)2 

'  AH) 

To  evaluate  the  interaction  matrix  elements,  we  use 
the  following  relation: 


1 

r 


/OO 

el'?(Ci-C2)K0(R|2gsin[(iy31  -  y>2)/2]| )dq. 

-OO 


(A4) 


Here  we  derive  the  rate  equations  and  the  stationary 
current  corresponding  the  the  setup  depicted  in  Fig[l0| 
a).  At  zero  magnetic  field  and  large  detuning  the  sin¬ 
gle  particle  ground  state  consists  of  the  Kramer  doublet 
\L1K'  f)  and  \L1K  J,)  localized  in  the  left  dot.  The  two- 
particle  states  are  characterized  by  their  charge  degree  of 
freedom  and  by  their  spin  valley  symmetry.  The  |2L,  OR) 
and  \1L,  1R)  states  with  antisymmetric  spin- valley  part 
are  non-degenarate  and  are  called  |  M)  and  |  M)  respec¬ 
tively.  They  are  defined  in  Eq.  ©  and  (fl5]) .  According 
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to  Fig{3  the  \2L,0R)  and  |1  L,  1 R)  states  with  symmetric 
spin-valley  part  are  each  three  fold  degenerate  and  in¬ 
clude  the  states  | FF)  and  IFF1)  respectively.  Following 
Eqs([T7])  and  ([23])  the  structure  of  these  states  is  given  by: 


with  Pi  +  Pp  +  Pp  +  Pm  P  Pm  =  1-  Elere  T e  denotes  the 
coupling  strength  to  the  emitter  contact  (right  contact 
in  FigjW]  a)  and  T c  denotes  the  coupling  strength  to  the 
collector/left  contact. 


| F_)  =  | FF)  =  |2 L,  0R}~  | K  |;  K  |) 
|F0)  =  \2L,0Rr}K  UK'  V+ 

I F+)  =  \2L,  0R)~  \K'  f;  K' t) 

I  F~)  =  \F~F)  =  \1L,1R)-\KUKI) 

|F0)  =  \lL,lRr  \K  UK'  V+ 

I F+)  =  \1L,1R}~\K'  t;/v'  f) 


Due  to  the  serial  setup  each  contact  couples  only  to 
its  adjacent  dot.  The  tunneling  rates  to  the  collector 
(left  contact  in  FigjTo]  a))  depend  on  the  correlated  two 
particle  states,  which  is  accounted  for  by  introducing  the 
spectral  weights  Sp  and  Sm  that  determine  how  easy 
an  electron  tunneling  between  the  left  dot  and  the  left 
reservoir  can  cause  a  transition  between  a  two  particle 
and  one  particle  state  on  the  dot. 

Sf  =  ^2  |  (L1K  j.|  ClLmKl  |F_)  |2 

m 

Sm  =  |  (L1K  ||  Q>LmK'\  | M)  | 2 


Due  to  the  simple  form  of  the  |1L,  li?)  wavefunctions 
the  spectral  weights  corresponding  to  the  tunneling  be¬ 
tween  right  contact  and  right  dot  has  trivial  spectral 
weights  Sm  =  \  (L1K  i\  amK^  \M)  \2  =  §  and  SF  = 
\{L\K  l\aR1Ki\F-)\2  =  1. 

We  denote  the  reduced  density  matrix  for  the  dou¬ 
ble  dot  by  P  and  denote  off-diagonal  matrix  elements 
with  P®  :=  (s'|P|s)  =  (Pss  )*,  where  s  and  s'  are  eigen¬ 
states  of  the  double  dot  and  diagonal  elements  with 
Ps  :=  (s|P|s).  We  neglect  processes  that  can  mix  states 
with  different  spin- valley  parfP^S  since  hyperfine  inter¬ 
action  is  small  in  carbon  nanotube  and  spin-orbit  cou¬ 
pling  is  mainly  a  spin-valley  coupling  with  a  very  weak 
effect  on  the  longitudinal  wave  function.  Then  degen¬ 
erate  states  are  equally  populated  and  the  relevant  ma¬ 
trix  elements  of  the  reduced  density  matrix  for  the  dou¬ 
ble  dot  system  are  given  by  Pi  =  P\lik'\)  +  P\lik\)\ 

PF  =  P\F+)+P\F0)  +P\F _>;  PF  =  PfI  +  PF°  +PfI  >  etc- 
The  corresponding  rate  equations  are  given  by: 

dtPi  =  -2TeP1+YcSfPf  +  2TcSmPm  (Bl) 

dtPF  =  —TcSfPf  —  itS(P£  -  pf) 

dtPp  =  Z-TEPi  +  its{Pf-P%) 

dfPp  =  (~-YcSf  +ie)PE  —  its{pE  —  pF) 

dtPM  =  -2T  cSMPM-itAs(pti-pM) 

dtPM  =  \vEPi+itAS{PV  -  P™) 

dtP]\i  =  — r  cSmPf  —  Has(Pm  —  Pm) 


These  rate  equations  are  solved  analytically.  The  re¬ 
sult  for  the  stationary  current  and  the  occupation  of  the 
states  are  given  by: 

Ion  =  =  ^IGSfSm^asTc/N  (B2) 

Ppp  =  3SMt2ASTE(4t2s  +  T2cS2F  +  4e2)/N 

N  =  12S m^as^ P  P  3SpSMt\sp%?pE 
+2SEtg(4:SMtAsPc  +  2t\sT  E  +  S\iT'E;TE) 

Here  tAs  denotes  the  interdot  tunneling  between  states 
|  M)  and  |  M)  and  ts  the  interdot  tunneling  between  the 
the  states  | FF)  and  | FF).  Spectral  weights  Sp  and  Sm 
account  for  the  overlap  between  an  electron  hopping  onto 
the  dot  and  the  two  particle  eigenstates,  and  Tc  and  T E 
are  the  coupling  strength  to  emitter  and  collector. 


The  stationary  current  (B2)  vanishes  as  1/e2  for  large 
e.  We  note  that  for  reversed  transport  voltage  still  a 
finite  current  can  pass  through  the  double  dot,  so  that 
for  large  e  the  double  dot  acts  as  a  current  rectifier.  We 
calculate  the  current  for  the  unblocked  voltage  direction. 
Inverting  the  bias  voltage  depicted  in  Fig  10  a)  current 


flows  by  electrons  hopping  from  the  left  reservoir  through 
the  double  dot  to  the  right  reservoir.  Applying  again  the 
rate  equation  approach  and  obtain: 


r  ^Rp  = _ 2eT%t\s _ 

°ff  h  1  h2t\s{TR  +  2SmTl)  p  smtlyr  1 


We  note  that  corresponds  to  Fc  in  Eq.(B2)  and 
corresponds  to  T ^ 
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